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This work is a continuation of our recent series of papers on Casimir 
friction, for a pair of particles of low relative particle velocity. Each 
particle is modeled as a simple harmonic oscillator. Our basic method, 
as before, is the use of quantum mechanical statistical mechanics, in- 
volving the Kubo formula, at finite temperature. In this work we begin 
by analyzing the Casimir friction between two particles polarizable in 
all spatial directions, this being a generalization of our study in EPL 
91, 60003 (2010), which was restricted to a pair of particles with lon- 
gitudinal polarization only. For simplicity the particles are taken to 
interact via the electrostatic dipole-dipole interaction. Thereafter, we 
consider the Casimir friction between one particle and a dielectric half- 
space, and also the friction between two dielectric half-spaces. Finally, 
we consider general polarizabilities (beyond the simple one-oscillator 
form), and show how friction occurs at finite temperature when finite 
frequency regions of the imaginary parts of polarizabilities overlap. 
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1 Introduction 



Casimir friction - a subclass of the Casimir field of research - has emerged to 
be a topic of considerable current interest. It is basically a non-contact kind 
of friction that may be related to electromagnetic fluctuations. The effect is 
small under normal circumstances, and may moreover be difficult to deal with 
theoretically since it involves energy dissipation, necessitating in turn the use 
of complex valued permittivities in macroscopic electrodynamics. Once one 
leaves the state of thermal equilibrium one will have to face fundamental 
problems. Thus, in connection with macroscopic electrodynamics a very 
instructive way of approach is to make use of the spectral summation method, 
but this method is based upon real eigenfrequencies in the spectral problem 
and how to deal with such a system in the case of dissipation is at present 
unclear. Cf., for instance, the discussion on this point in Refs.JI] and [2]. 

Most of the previous works in this area are based on the macroscopic 
dielectric model with permittivity properties. See, for instance, Refs. [21 HJ 
El El El El E] • There exists, however, a different strategy in order to deal with 
Casimir friction, namely to consider the statistical mechanics for harmonic 
oscillators at finite temperature T moving with constant velocity v relative 
to each other. We argue that such a microscopic model, in spite of its sim- 
plicity compared with a full-fledged macroscopic model, has nevertheless the 
capacity of providing physical insight in the problem. And this is the kind 
of model that we shall consider in the following. We shall make an exten- 
sion of our considerations in recent papers [IH1 [HI [12]; cf. also the arXiv 
report [13]. And this is again based upon a study of ours back in 1992 [14J, 
dealing with the same kind of system. It contains a generalization to the 
time-dependent case of the statistical mechanical Kubo formalism spelled 
out for the time- independent case in Ref. [15] . 

Whereas our previous works were limited to the case of two microscopic 
oscillators, we shall here generalize the formalism so as to deal with one 
oscillator outside a slab (a collection of oscillators). Also, the generalization 
to two slabs in relative constant motion is easily achievable. 

We emphasize that we are considering dilute media only. The general- 
ization to media of arbitrary density ought to be tractable with the actual 
method, but has to our knowledge not been treated yet. Our results are 
that the energy change AE because of Casimir friction is finite in general. 
This corresponds to a finite friction force. At zero temperature the formal- 
ism yields AE — > 0, however; this being due to our assumption of a slowly 
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varying coupling. For rapidly varying couplings, there will be a finite friction 
force also at T = PUDS]- 

We mention finally that the microscopic approach has been analyzed by 
other investigators also, especially by Barton in recent papers [HI [T71 IT8] . 
The equivalence between Barton's results and our own results is not so easy 
to see by mere inspection since the methods are different, but were shown 
explicitly to be equivalent in one of our recent papers [T2] . 



2 Friction force between two oscillators 

Consider the quantum mechanical system of two polarizable particles whose 
reference state is that of uncoupled motion corresponding to a Hamiltonian 
H . The equilibrium situation is then perturbed by a time independent 
term —AF(t) where A is a time independent operator and F(t) is a classical 
time dependent function. For a pair of polarizable particles perturbed with 
dipole-dipole interaction we have 

- AF(t) = ipijS U s 2 j, (1) 

where the summation convention for repeated indices i and j is implied. Here 
s\i and s 2 j are the components of the fluctuation dipole moments of the two 
particles = 1,2,3). With the electrostatic dipole-dipole interaction, 

*» = - - . P) 

where r = r(t) with components Xi = Xi(t) is the separation between the par- 
ticles. In an earlier work we also studied the situation with time-dependent 
interaction (with retardation effects) [Tj5], but we will avoid this added com- 
plexity here. 

The situation with which we shall mainly be concerned in the follow- 
ing is when the relative velocity v of the oscillators is constant. Then the 
interaction will vary as 



AF(t) 



Sli$2j, (3) 



and the components of the force B between the oscillators are 

d 
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The equilibrium situation with both particles at rest is represented by the 
first term in (J2J). It gives rise to the (reversible) equilibrium force. Thus the 
friction is connected with the second term. To simplify, we shall here neglect 
the first term by which the two oscillators will be fully uncorrelated in their in 
their relative position r = Tq. Thus we can write — AF(t) — > —AiFi(t) where 
Ai = B\ and Fi(t) = V[t. The friction force will be a small perturbation upon 
the equilibrium situation, and it leads to a response A(Bi(t)) in the thermal 
average of B\. According to Kubo [20j [151 E] 

/oo 
<t>BAl q {t-t')F q {t')dt\ (5) 
-oo 

where the response function is 

<j>BAi q {t) = ^TrpL4 9 , (6) 

Here p is the density matrix and Bi(t) is the Heisenberg operator Bi(t) = 
e itH/h Bi e - itH l h where like A is time independent. With Eqs. © and (0} 
and with Fi(t) = V[t expression (6) can be rewritten as 

4>BAlq(t) = Gl g ij nm (j)ij nm (t), (7) 

where 

__ d%_ (9Vw 
lqijnm dx x dx q ' 1 ' 

4>ijnm(t) = Tl{pC ijnm (t)}, (9) 
Cij nm (t) = — [SuS2j, Si n (t)s2m(t)] (10) 

(the i in the denominator is the imaginary unit). 

Here as in Refs. [H] and [10] the perturbing interaction ([Q) will be con- 
sidered weak. This will also hold for dilute dielectric media by which thermal 
averages of products, containing su and s 2j as factors, to leading order fac- 
torize. Further, with isotropy or scalar polarizability the components of each 
of the fluctuating dipole moments are uncorrelated. 

From now on we find it convenient to utilize imaginary time, which was 
used in Sec. 4 of Ref. [TJ]. With this we get 



9ijnmW = T^[pSln(t)s 2m (t)s U S2j] = 5'lm(A)fi'2jm(A), (11) 
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9apqW = (s aq (t)s ap ) = g a (X)S qp , (a = 1,2), (12) 

where g(X) is the correlation function and angular brackets denote thermal 
averages ((..) = Tr[p...]). The A is imaginary time given by 

\ = i t -, (13) 

so that for an operator B 

B(t) = e XH Be- XH . (14) 

With this, 

4>ijnm{t) = TjT [9ijnm (P + A) - (A)] 

= ^ W + A)^ 2 (/3 + A) - gi{X)g 2 (X)} 5 in 5 jm . (15) 

The response function 4>ij nm (t) corresponds to the retarded Green function, 
in the usual language of quantum field theory 

Earlier we unfortunately made a mistake by defining A in Eq. (fTB"]) with 
opposite sign [14] . However, this did not influence the results of previous 
applications since the operators A and B were equal apart from prefactors. 

In Appendix B of Ref. [H] it was found that 

4>(co)=g(K), (16) 

so that the Fourier transforms of the response function <fi and the correlation 
function g are equal. Explicitly, 

poo 

4>(cj) = / ct){t)e- iult dt ((pit) = for t < 0), 
Jo 

~g(K)= g(X)e iKX d\, (17) 

and 

K = ihw. (18) 

[ In Ref. [H] K = —ihu was used due to the mistake mentioned.] Equality 
( FTBl holds in the common region |Im(i^)| < C with C > and Im(tu) < 0, 
i.e., Ke(K) > where both functions are analytic. 
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It may be remarked that Eqs. (1121) - (1181) above correspond in a quantum 
field theoretical language to the statement that the spectral correlation func- 
tion (frequency u) is equal to the imaginary part of the spectral retarded 
Green function multiplied with coth(^(3hu). See, for instance, Eq. (76.6) 
in Ref. [21J. In the present case the function g+(ou) which is the real time 
Fourier transform of g+(t) = ^(g(it/h) + g(—it/h)), will be the spectral cor- 
relation function while the <p(—uj) (with Fourier transform f lTTj) ) will be the 
spectral retarded Green function (for u > as g+(ou) is symmetric in u while 
lm[0(— u)] is antisymmetric). 

With Eqs. (HU) and (d2D we have 

(A) = g(X)S in Sj 

g(X) = gi (X)g 2 (X). (19) 
Thus in K- space the g{K) can be written as the convolution 

9( K ) = ^Y,9i( K o)h(K-K ), (20) 

which is Eq. (4.15) of Ref. [H]. 

An advantage of using imaginary time is that g a {K) can be identified 
with the frequency dependent polarizability a a K of oscillator a (= 1,2). For 
a simple harmonic oscillator with eigenfrequency u a one has (15] 



9 ° (A) = a - = 7^T(^- (21) 

where a a is the zero-frequency polarizability. [ The gi(K) given by Eq. (4.10) 
of Ref. [E] will differ from the one of Eq. ( J2TT) by a factor e 2 where e is the 
electron charge since here the s a is identified with a dipole moment.] 
With the above one finds for the Fourier transform of ([9]) 



[u) = <f>(u))S iri 5 jm , (22) 



with <f)(u) = g(K). This follows by use of Eqs. f|T5l) and f|T9l) from which 
gijnm(K) = g(K)5 in Sj m , and when account is taken of Eq. f|T6l) . The g(K) 
is moreover given by the convolution (120]) where for two simple harmonic 
oscillators the g a (K) (a = 1, 2) are given by Eq. (|2~T|) . 
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To obtain the perturbing force (EJ) the expression ([7]) should be evaluated. 
With (12 2D one can write 



4>ijnm{t) = 4>(t)S in 5 



j mi 



by which 



4>BAlq{t) = Gl q (j)(t), 
Glq Glqijij TnjT q ij } 

dip i:j 



T, 



lij 



dxi 



With ipij given by Eq. (|2]) one finds 



T, 



lij 



YhXiXjXi 3(xi5ij + XjS u + xAj) 



(23) 

(24) 
(25) 

(26) 
(27) 



by which 



Glq TnjTqij 



+ 



270 

r 12 



r + 



X[Xq 



, I ^5lq = ^5lq + 



36xix 



iXq 



,10 



(28) 



With the above expressions the result for the perturbing force ([5]) will be 
precisely the same as obtained in Ref. [10], except from subscripts / and q. 
One may simply insert Eq. fl24l into Eq. ([5]) with Gi q given by Eq. (128]) . 4>(t) 
following from the relations (fT6]) - (l23l . 

Thus, like in Eqs. (8) and (9) of Ref. [TO] we here get the perturbing force 



where 



F l = A(Bi(t))=F r i + F fh 

GlgVqt j (j)(u)du 



r! 



(29) 



(30) 



is part of the reversible force. The part of the force representing friction is 

POO 

Ffi = -GiqVq / 4>{u)udu. (31) 

JO 
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The Fourier transformed version of this equation is, like Eq. (11) of Ref. [TO]). 

d(j>(u)) 



Ffl = -lGl q V q - 



uj=0 



(32) 



The expression for </>(u;) is evaluated below Eq. ( 1581 in the next section, by 
which the explicit expression (14*51) is obtained for Ff\. 

Now the dissipated energy may be obtained; this requiring a perturbing 
interaction to last a finite amount of time to make it unique. Thus the force 
expression may be modified to 

Fi(t) -> qi (t), (33) 

where qi(t) is interpretable as a position. The dissipated energy AE d is then 
given by expression (27) of Ref. [TU|. In the present case this will be (Ai = B{) 



/oo |~ rt 

qi(t) / <f>AAi q {t - t')q q {t')dt' 
-oo J —oo 



t 

' dt. (34) 



3 Friction between particle and a half-space, 
and between two half-spaces 

Armed with the result (l2"8l for Gi q it is now straightforward to extend the 
results in the previous section to the situation where a polarizable particle 
moves parallel to a resting dielectric half-space. We assume then that the 
particle density p in the half-space is low so that the forces are additive. Let 
the half-space be located at z > z such that its surface is parallel to the 
xy plane at the vertical position z = Zq. The dielectric particle is located at 
the origin and moves with constant velocity v x = v\ along the x axis. The 
resulting friction then follows by integrating Gi q with I — q — 1 over the 
dielectric half-space. With expression (128]) we find 



Gh = p Gndxdydz. (35) 

J z>zo 

Symmetry with respect to the x and y coordinates means that the x 2 in 
Eq. (1281) can be replaced with ^(x 2 + y 2 ), and we can use cylindrical coordi- 
nates with p 2 = x 2 + y 2 and dxdy = l-npdp. Thus with r 2 = p 2 + z 2 , 

Gxxdxdy = 36tt j + pdp = 0, ( 36 ) 
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by which 



CO 



157T 3np 



G i = Pl ^dz = ^. (37) 



o~6 n 7 5 

z zz zz 



Finally, for two dielectric half-spaces moving parallel relative to each other 
one can obtain the friction force per unit area from 

/°° 3n 
G h dz = ^jPip2, (38) 

where d is the gap width and p\ and p 2 are the (low) particle densities in the 
two half spaces. 

The function <p(u) = g{K) is needed to obtain the friction. Performing the 
summation in Eq. (1201) . inserting the expression (12 ip we obtain, in agreement 
with Eq. (4.16) in [H], 

~g(K) = Hf(K), (39) 

with 

E 1 E 2 a 1 a 2 

ti 



4smh(±/3£ 1 )sinh(±3£ 2 )' 
, s Eisinh(i/3Ei) £ 2 sinh(±/3£ 2 ) , s 

™= wrn + v + % - (40) 

Here Si and £ 2 are defined as 

Ex = + £ 2 , E 2 = £ x - £? 2 , = fiwi (i = 1, 2). 

When the velocity is small and constant (or very slowly varying), only the 
limit K — Y is needed. One further sees that the contribution requires that 
E 2 — > 0. Because of this, the f(K) becomes a 5-function (plus a constant) 

f{K) = ~PK6(Z 2 ), (41) 

like Eq. (4.18) of Ref. [H]. [Here Re(f^) > in view of the correction 
mentioned below Eq. ffl5|) .] To obtain the friction force like ( 132]) for a pair 
of particles and thus the more general situation, the derivative of f(K) with 
respect to u is needed. With <5(E 2 ) = S(H(cji — w 2 )) = 5(u>i — io%)jh and 
K = ifvjj one finds 

^ /W = _ S ^) = | Wwi _ W2 ). (42) 
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For a pair of polarizable particles the friction force (I3"2"j) thus becomes 



F fl = -G lq v q H^-5{u x -uj 2 ), (43) 



with H given by Eq. (j4"0|) and Gi q by Eq. 

For the more simple model studied in Ref. [10] the Gi q v q is replaced by 
G = (V"?/>)(v • Vip). The result (19) of that reference is recovered if the 
{tvjj a ) 2 a a (a = 1,2) is replaced with h 2 /m a in the expression (jJOJ) for H. 
This replacement follows from a corresponding change in g a (K) as given by 
Eq. ([UJ) to the one given by Eq. (4.10) in Ref. [H]. 

It is now straightforward to obtain the friction between a polarizable 
particle and a half-space. One can simply replace the G[ q v q in Eq. fj43|) with 
GhV , where t> is the velocity of the particle parallel to the plane and Gh is 
given by expression f )3T|) . Thus 

F. =-^1^(^-0*). (44) 
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Likewise, with two half-spaces moving relative to each other the G/j is re- 
placed with G given by Eq. (|38|) to obtain the friction force per unit area 

^ 7T -pip 2 H^-5(uj 1 - u 2 ). (45) 



4 General polarizability 

For simple harmonic oscillators the polarizability is given by Eq. (12 ip . How- 
ever, it can be a more general function of K that may be regarded as resulting 
from a sum of harmonic oscillators. Thus we may write 

h(K 2 ) = UK) = a aK , (46) 

where it can be shown [22J that the function h(K 2 ) satisfies the relation 

M^) = / ^^(>* 2 ), («) 

with 

a a (m 2 )m 2 = lm[h(— m 2 + ij)), (m = fku = —iK, 7 — > + ). (48) 
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This is analogous to the Kramers-Kronig relations for the dielectric permit- 
tivity. An implication of the more general expression for g a (K) is that the 
a a (a = 1,2) in Eq. (J1D]) is replaced with J a a (m^)d(m^). Integration of the 
5-f unction in Eqs. fj42|) - P5|) further gives 



J 6(u}t - u 2 )d(m 2 1 )d(ml) = 4(^w 1 ) 2 'duj^ (49) 

So altogether the H(ti(3/2)5(ui — oj 2 ) in these equations will be replaced by 

(u>i, (jJ 2 — >■ ui,m = fvjj) 



H Q = -£- I (5°) 



7r/3 r m A ai(m 2 )a2(m 2 ) 
T J sinh 2 (i/3m) 

With this the friction forces, Eqs. ([4"5]) -(l4"5 ]) . for a pair of particles, one particle 
and a half-space, and two half-spaces, become respectively 

Ffi = —Gl q V q H , 



3n p 

pip2H - (51) 



Fh — tH() 

3tt 



8d 4 ' 

Thus, at a finite temperature a finite friction force is obtained when finite 
values of «i(m 2 ) and a 2 (m 2 ), as defined in Eqs. (HB|) - (H5|) . overlap. 
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